Crossover temperature of Bose-Einstein condensation in an atomic Fermi gas 
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We show that in an atomic Fermi gas near a Feshbach resonance the crossover between a Bose- 
Einstein condensate of diatomic molecules and a Bose-Einstein condensate of Cooper pairs occurs 
at positive detuning, i.e., when the molecular energy level lies in the two-atom continuum. We 
determine the crossover temperature as a function of the applied magnetic field and find excellent 
agreement with the experiment of Regal et al. [Phys. Rev. Lett. 92, 040403 (2004)] that has recently 
observed this crossover temperature. 



Introduction. — An atomic Fermi gas near a Fesh- 
bach resonance is a fundamentally new superfluid sys- 
tem. The reason is that near a Feshbach resonance the 
gas does not only consist of atoms but also of diatomic 
molecules. Moreover, the energy difference betvifeen the 
molecular level and the threshold of the two-atom con- 
tinuum, known as the detuning 5, can be experimentally 
tuned by means of a magnetic field In combination 

with the fact that for fermionic atoms these molecules 
are very long lived H, ^ H, 111 ; such a gas thus offers the 
exciting opportunity to study in detail the crossover be- 
tween the Bose-Einstein condensation (BEG) of diatomic 
molecules and the Bose-Einstein condensation of atomic 
Cooper pairs, i.e^ the Bardeen-Cooper-Schrieffer (BCS) 
transition [3, S S ll^ ■ Indeed, a Bose-Einstein conden- 
sate of molecules has recently been observed |0 0, • 
More recently, a claim for Bose-Einstein condensation 
of atomic Cooper pairs was made We will show, 

however, that the date reported in Ref. Q can be 
understood in terms of a Bose-Einstein condensation of 
molecules. 

At zero temperature the physics of the BEC-BCS 
crossover occurring near a Feshbach resonance can be 
understood as follows. The superfluid phase of the gas 
is always associated with a Bose-Einstein condensate of 
pairs of atoms, but the wave function of the pairs is given 
by 



v/z(^Xm(x,x') I closed) + VT^2X^Xaa(x,x';5)|open) . 

At large negative detuning the energy of the molecule lies 
far below the threshold of the two-atom continuum and 
we have Z{5) ~ 1. In that case we are dealing with a 
Bose-Einstein condensate of diatomic molecules and the 
spatial part of the pair wave function is equal to the 
(bare) molecular wave function Xm(x, x'). The spin part 
of the pair wave function is then equal to |closed), i.e., the 
spin state of the closed channel of the relevant Feshbach 
problem At large positive detuning the molecular 

energy level lies far above the threshold of the two-atom 
continuum and can be (adiabatically) eliminated. We 
then have that Z((5) ~ and the spatial part of the pair 
wave function equals the usual BCS wave function for 
atomic Cooper pairs Xaa(x, x'; 5). This Cooper-pair wave 



function depends on the detuning, because the effective 
attraction between the atoms depends on the detuning. 
The spin state of the Cooper pairs is, however, always 
equal to the spin state of the open channel of the Fesh- 
bach problem, denoted here by |open). 
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FIG. 1: Phase diagram of atomic ''"K as a function of mag- 
netic field and temperature for a Fermi temperature of the 
gas of Tf = 0.35 /iK. The solid line gives the critical tem- 
perature for either a Bose-Einstein condensation of molecules 
or a Bose-Einstein condensation of atomic Cooper pairs. The 
critical temperature for the latter is calculated by simultane- 
ously solving the BCS gap equation and the equation of state 
of an ideal mixture of atoms and molecules. In this equation of 
state, we use the energy of the molecules given by Eq. I^J be- 
low. For comparison, the upper dashed curve is the analytical 
BCS result T /Tp =j;8e^"V'r)e-''''^''Fl''l , where 7 ~ 0.5772 is 
Euler's constant |7l[lq|. The lower dashed line is the crossover 
between the two Bose-Einstein condensed phases, which is the 
main topic of this Letter. 

With this physical picture in mind, the point where 
the crossover takes place is thus determined by the de- 
tuning at which the amplitude Z{5) relatively abruptly 
crosses over from one to zero. Based on two-body physics 
[isf we would expect the crossover to occur exactly on 
resonance, i.e., at zero detuning. It is the main purpose 
of this Letter, however, to point out that many-body 
physics changes this result and shifts the crossover point 
to positive detuning. Quantitatively, the crossover oc- 
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curs, for an incoherent mixture of two hyperfine states 
with equal density n/2, at the detuning where the molec- 
ular energy level becomes equal to twice the Fermi energy 
ep — ft^(37r^n)^/'^/2m. Using the theory presented below 
we can also accurately determine the crossover point at 
nonzero temperatures. The resulting (mean-field) phase 
diagram for atomic '^''K under the conditions of the exper- 
iment of Regal et al. is shown in Fig. ^ and summarizes 
the main conclusion of our work. 

Note that in the usual BEC-BCS crossover problem, 
studied in condensed-matter physics in the context of 
the high-temperature superconductors, Z{S) is always 
identically zero and the crossover is associated with a 
qualitatively different behavior of the Cooper-pair wave 
function Xaa(x, x'; 5) This emphasizes the 

fundamentally new nature of the superfluid state in an 
atomic Fermi gas near a Feshbach resonance, which actu- 
ally shows a macroscopic coherence between atoms and 
molecules. 

Poor man's approach. — Before discussing the the- 
ory that incorporates the resonant interactions between 
atoms, we first consider the case of an ideal mixture of 
molecules and atoms to establish the physical picture of 
the crossover most clearly. At positive detuning a sta- 
ble molecular state does not exist, because the molecule 
can energetically decay into two free atoms. In previ- 
ous work |l5l Il9|. however, we have shown that in first 
instance it is reasonably accurate to neglect the finite 
lifetime of the molecule. Close to the Feshbach reso- 
nance, the interaction with the atomic continuum shifts 
the molecular energy level downward from the detuning 
5to em =i h'^/ma^. Here a{B) = ai,g[l - AB / {B - Bq)] is 
the full atomic s-wave scattering length of the Feshbach 
resonance, which is experimentally characterized by its 
location at magnetic field Bq, its magnetic field width 
AB, and the background scattering length abg. 

If the molecules are Bose-Einstein condensed, the 
chemical potential of the atoms is equal to em/2 — 
h^/2ma^. If wc take the atoms to be noninteracting, 
we can thus at zero temperature immediately determine 
the density of atoms in the Fermi sea below the chemical 
potential. Subtracting this result from the total density 
n, and dividing by two, we find that at zero temperature 
the density of molecules in the Bose-Einstein condensate 
is equal to 
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(fcF|a|)= 



(1) 



where kp = {Sn^nY^^ is the Fermi momentum of the gas. 
Note that this estimate is only valid for positive detuning, 
where a < 0. So at resonance the density of condensed 
molecules is just n/2, i.e., half the total density of atoms 
in the gas. Moreover, the density of condensed molecules 
vanishes at kp\a\ — 1. Physically, this situation occurs 
when the energy level of the molecule is exactly equal to 
twice the Fermi energy of the gas. This result is sensible. 



since if the energy level of the molecule is higher, there 
will be no molecules at zero temperature. The whole 
gas then consists of atoms. In the experiment of Regal 
et ai, this condition gives a magnetic field of 0.5 Gauss 
above the resonance, in excellent agreement with the data 
shown in their Fig. 2. 

To find the same criterion at nonzero temperatures 
is also possible. We know that for temperatures below 
the Fermi temperature the density of fermionic atoms is 
hardly influenced by temperature. Not too close to reso- 
nance the molecular condensate density is just depleted 
by thermal fluctuations, i.e., a thermal cloud of molecules 
forms with increasing temperature. Calculating the crit- 
ical temperature for a density of ideal Bose molecules 
given by Eq. Q gives us the result 
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with Tp = ep/fce the Fermi temperature and ^(3/2) ~ 
2.612. This result can be directly compared with the data 
of Regal et al. presented in the x-y plane of their Fig. 4. 
In view of the simplicity of the approach, the agreement 
is remarkable. 

Molecular Bose-Einstein condensate. — To properly 
incorporate the resonant interactions between the atoms, 
a more involved treatment of the gas is necessary. In- 
troducing creation and annihilation operators for the 
molecules and atorns, the grand-canonical hamiltonian 
of the gas becomes |20l | 



H 



4m 



2/i 



+ J2 J dx^pli^) 



2m 



(3) 



dxg ^ipl,{x)'ip'^{x)'ipi{x) -)- V|(a;)V'j(a;)V'm(a;) 



where the two hyperfine state of the atoms are de- 
noted by I t) and | J,), the atom-molecule coupling 
constant g = hy/ inahgABAfi^^g^g/m, and the magnetic 
moment difference A/imag between the hyperfine states 
|open) = (I tl) — I J,t))/'s/2 and |closed) gives the de- 
tuning 6 — A/imag (B — Bo) Finally, the molecular 
energy is approximated by the energy where the molec- 
ular density of states has a maximum. To find this max- 
imum we use that in the two-body T-matrix approxima- 
tion the frequency-dependent self-energy of the molecules 
is —irjVhio [lajllyl- I'^ general, this gives 



(4) 



where 77^ — g'^m^ /Wir^hP is the energy scale associated 
with the width of the Feshbach resonance. This energy 
scale is in fact of fundamental importance, because it 
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shows that at zero temperature the thermodynamic prop- 
erties of a resonant atomic Fermi are not solely deter- 
mined by the Fermi energy. This is particularly true for 
the experiment of Regal et al., for which 77^ ^ ep, since 
for the Feshbach resonance of interest rj^ ~ 7.7 mK. Close 
to resonance, where 5 -^rf , the molecular energy reduces 
to /rf^ , which can be shown to be equivalent to Y? jma^ 
as expected. Note that the omission of the background 
scattering length in the hamiltonian is justified because 
the region of interest takes place relatively close to reso- 
nance. 

To find the crossover temperature for positive detun- 
ing, we consider the gas to have a Bose-Einstein conden- 
sate of molecules, and perform a quadratic expansion of 
the hamiltonian around the nonzero expectation value 
(^m(x)) = ^rimc- This leads to the ideal gas expression 
for the molecular density 

"■m = "mc + y X! g£k/2fcBT _ ^ ' 
k 

where V is the volume of the gas and ek = 'f?\<? j2m. 
However, for atoms with momentum k, the resulting 
hamiltonian leads to a fluctuation matrix 

ek - em/2 5V"mc 
ffV^W -(ek-ein/2)J ' 

which can easily be diagonalized by means of a Bogoli- 
ubov transformation. Performing the calculation, we ul- 
timately find for the total atomic density 




where the dispersion for the atoms obeys = 
\/ (ek - em/2)2 -I- 5^nmc- 

For fixed positive detuning and temperature, the equa- 
tion for the total atomic density n — In^-a+n^ determines 
the molecular condensate nmc- The result of these cal- 
culations for the experiment of Regal et al. is shown in 
Fig. 121 This figure can be directly compared with the 
data in their Fig. 4. Again the agreement is remark- 
able. Having said that, it is important to realize that 
the experiment is performed in an optical trap, whereas 
we have considered the homogeneous situation. General- 
izing our poor man's approach to the trapped situation 
shows, however, that this does not affect the position 
of the crossover line, because in that case the homoge- 
neous criterion is satisfied in the center of the trap. The 
inhomogeneous analysis can be carried out in the local- 
density (or Thomas- Fermi) approximation, but this is be- 
yond the scope of the present paper and is left for future 
work. Such an analysis is certainly needed to obtain a full 
quantitative agreement between theory and experiment. 




FIG. 2: Molecular condensate fraction rime /rim in an atomic 
function of magnetic field and temperature for 
a Fermi temperature of the gas of Tf = 0.35fiK. This figure 
should be compared with Fig. 4 of Ref. p^. 



The most important approximation that we have made 
in our calculation of the crossover temperature is to ne- 
glect the finite lifetime of the molecules. Including this 
finite lifetime is not an easy task, because a consistent 
approach requires that the self-energy of the molecules 
is calculated at least in the many-body T-matrix ap- 
proximation, just as we have done in our recent work 
on the observation of molecular Kondo resonances in an 
atomic Fermi gas near a Feshbach resonance The 
physical reason for this complication is that the decay 
of the molecules can be Pauli blocked by the presence 
of the atomic Fermi sea. It is this Pauli blocking that 
is ultimately responsible for the molecular Kondo reso- 
nances and it will, therefore, also play an important role 
in a quantitative analysis of the molecular lifetime effects. 
Qualitatively, we expect that a broadening of the molec- 
ular energy level will not have a substantial effect on the 
location of the crossover line. It will, however, lead to an 
increase of the molecular condensate fraction, because of 
the presence of molecular states with energies below Cm- 
We believe that these lifetime effects may also have a 
bearing on the considerable narrowing of the molecular 
thermal cloud that was also observed by Regal et al. |Q| . 
Work in this direction is presently being completed and 
will be reported elsewhere. 
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